A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 

08.16.3; 08.01.1; 08.01.3; 08.09.3; 02.12.3 



ASTRONOMY 

AND 
ASTROPHYSICS 



The beryllium abundance in the very metal-poor halo star 
G 64-12 from VLT/UVES observations* 



o- 
o: 
o. 

(N ■ 

CW 

<D ■ 

C/V 

OV 

<n : 



F. Primas 1 **, M. Asplund 2 , P. E. Nissen 3 , and V. Hill 1 

1 European Southern Observatory, Karl-Schwarzschild Str. 2, D-85748 Garching b. Munchen, fprimas@eso.org, vhill@eso.org 

2 Uppsala Astronomical Observatory, Box 515, SE-751 20 Uppsala, Sweden, martin@astro.uu.se 

3 Institute of Physics and Astronomy, Aarhus University, Denmark, pen@ifa.au.dk 

Received ; Accepted 



> 

00 . 

ov 
o; 
o. 

Or 
i ' 

O' 

a: 

c/a . 



* 



Abstract. We report on a new spectroscopic analysis 
of the very metal deficient star G 64-12 ([Fe/H]=— 3.3), 
aimed at determining, for the first time, its beryllium con- 
tent. The spectra were observed during the Science Ver- 
ification of UVES, the ESO VLT Ultraviolet and Visible 
Echelle Spectrograph. The high resolution (~48 000) and 
high S/N (~T30 per pixel) achieved at the wavelengths 
of the Be II resonance doublet allowed an accurate de- 
termination of its abundance: log N(Bc/H) = —13.10 ± 
0.15 dex. The Be abundance is significantly higher than 
expected from previous measurements of Be in stars of 
similar metallicity (3D and NLTE corrections acting to 
make a slightly higher value than an LTE analysis). When 
compared to iron, the high [Be/Fe] ratio thus found may 
suggest a flattening in the beryllium evolutionary trend at 
the lowest metallicity end or the presence of dispersion at 
early epochs of galactic evolution. 

Key words: Stars: Population II, abundances, atmo- 
spheres - - Galaxy: halo - - Instrumentation: spectro- 
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1. Introduction 

The knowledge of lithium, beryllium, and boron in stars 
play a major role in our understanding of Big Bang nu- 
cleosynthesis, cosmic-ray physics, and stellar interiors. 
Lithium (together with D and He) is a key element to 
probe the Big Bang nucleosynthesis scenario, because 
most of 7 Li is primordially produced. On the contrary, 
6 Li, 9 Be, 10 > 1:L B (with a possible contribution to n B from 
z/-spallation in supernovae being still under debate) were 
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suggested to originate primarily from spallation reactions 
in the interstellar medium (ISM) between cosmic-ray (CR) 
a-particles and protons and heavy nuclei like carbon, oxy- 
gen, and nitrogen (Reeves et al. 197C). The production 



rate of any of these three light elements is then propor- 
tional to the product of the abundance of CNO (propor- 
tional to the number of supernovae, N$n) and the CR 
flux (assumed to be proportional to the supernovae rate) . 
A slope of two in the logarithmic plane 6 Li, 9 Be, 10 B vs 
O is thus expected. If [O/Fe] in metal-poor halo stars is 
constant, then this same slope holds vs Fe as well. 
On these premises, the picture emerging from the first 
systematic analyses of Be and B in stars of the Galactic 
halo and disk from high resolution data clearly challenged 
this scenario: both elements scale almost linearly, and 
not quadratically with metallicity from the early Galactic 
ep ochs up to now (e.g. Molaroetal. 1997 , Boes gaard et 
al.~ 



1999b 



1997, Primas et al. 1999 for 



B 



for Be; Duncan et al. 
Such finding has been interpreted as an indication for 
a primary (instead of secondary) origin of Be and B, that 
in turn may imply the need for a production mechanism 
independent of the metallicity in the ISM. Although the 
most recent theoretical scenarios (which seem to converge 
towards the proximity of supernovae as one of the most 
likely and efficient sites for Be and B production) are able 
to reproduce the observed abundance levels and trends, 
new observational inputs are still in high demand in order 
to further constrain the models. The evolution/behavior 
of Be at the lowest metallicities, the presence and mag- 
nitude of scatter (if any), the presence and location of a 
possible change of slope along the trend, and a success- 
ful solution of the oxygen dilemma (high or flat) in the 
Galactic halo are some of the still missing inputs. 
With this Letter, we aim at adding a new tile to the field of 
light elements evolution, presenting the first beryllium de- 
tection in a star characterized by a metallicity well below 
one thousandth solar ([Fe/H]=— 3.3), the halo star G64- 
12. Our newly determined [Be/Fe] ratio will be compared 
to previous observational analyses in low metallicity stars 
and confronted with different theoretical suggestions. 
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Table 1. Observations log for G 64-12 



V Date Mode A c A-coverage Exp. Time S/N a 

nm nm min 

11.46 Feb. 12 Dic#l 346 (CD#1) 305-385 2x90 75,80 

580 (CD#3) 480-680 2x90 350,370 

Feb. 14 Dic#l 346 (CD#1) 305-385 2x80 65,70 

860 (CD#4) 665-1050 2x80 270,280 



measured on each single spectrum, around 313 nm, 510 nm and 700 nnr in the B346, R580, and R860 settings respectively 



2. Observations and Data Reduction 

A total of 6 hours integration time (cf Table 1 for more 
details) were devoted to G 64-12 (V=11.45 mag), during 
the Science Verification of UVES (February 10-18, 2000). 
The target was observed in Dichroic mode, and the six 
hours were splitted in four exposures of ~ 90 min each: 
two with the standard Dichroic #1 setting and two with a 
slightly modified Dichroic #1 setting. All the spectra were 
observed with a slit of 0.8" , corresponding to nominal re- 
solving powers of 48 000 in the blue and 55 000 in the red, 
later confirmed by the Thorium- Argon calibration frames. 
The spectra were reduced within the MIDAS UVES con- 
text. First, bias and inter-order background were sub- 
tracted from both science and flat-field frames, then the 
object was optimally extracted, divided by the flat-field 
(extracted with the same weighted profile as the star), A- 
calibrated and the extracted orders merged. All the spec- 
tra were then registered for barycentric radial velocity 
shifts, co-added, and normalized. An independent check 
was performed using standard IRAF tasks for echelle spec- 
tra reduction, which produced the same final product. The 
S/N ratio achieved at the Be II resonance doublet on the 
final spectrum (cf Fig. 1) is 130 per pixel. 

3. Stellar Parameters 



is too uncertain to be used to derive the surface gravity of 
the star. We note, however, that the position of G 64-12 
in the (b — y)o — cq diagram (see Ryan et al. 1999, Fig. 3) 



Using the Infrared Flux Method, Alonso et al. ( 1996] ) de- 
rived T cS = 6470 ± 90 K for G 64-12. This is consider- 
ably higher than the value of T e g= 6220 K determined 
by Ryan et al. (1999) from several colour indices using 
the "low" T cff scale of Magain fll987| ). It should be noted, 
however, that Ryan et al. adopt a negligible low interstel- 
lar reddening of G 64-12, whereas the Stromgren uvby-ji 
photometry as given in their Table 1 suggests a reddening 
E{b-y) = 0.023 corresponding to E(B - V) ~ 0.03. Cor- 
recting the colour indices for this reddening would increase 
the Ryan et al. temperature to about 6360 K. We conclude 
that the effective temperature of G 64-12 is quite uncer- 
tain; depending on the adopted reddening and the zero- 
point of the T e g scale, one derives values between 6200 
and 6500 K. However, fitting of the Ha wings indicates 
a temperature closer to the higher end of the considered 
range of values (T cff ~ 6400 K). 
The Hipparcos parallax of G 64-12 (w = 1.88 ± 2.90 mas) 



is similar to that of HD 84937. Using t he ab solute magni- 
tude calibration of Nissen & Schuster ( 1991 ), My — 3.6 is 
derived for both stars. Hence, the gravity of G 64-12 must 
be close to that of HD 84937 (log g= 4 .1) as derived from 
its Hipparcos parallax (Nissen et al. 1997 ). Altogether, 
we estimate that the gravity of G 64-12 lies in the range 
log g= 3.9 - 4.3. 

The metallicity of G 64-12 was derived from equivalent 
widths of 22 Fe I lines in the spectral region 480 - 560 nm. 
The lines are weak (2 < EW < 30 mA), but the EWs 
are very reliable due to the high S/N. Adopting ^/-values 
from O'Brian et al. ( 1991 ) and using a ID model atmo- 
sphere with T oS = 6400, [Fe/H] = -3.28 ± 0.08 dex is 
derived. Allowing for the possible range of T c g (6300 - 
6500 K), [Fe/H] = -3.28 ±0.10. 

It is worth noticing that the same procedure applied to 
5 Fell lines yields a somewhat lower iron content [Fe/H] 
— —3.42 ± 0.10. Though this would be in principle our 
preferred metallicity indicator (Fe n lines are much less 
sensitive to NLTE corrections), the gf- values of these lines 
are not very accurately known. Four out of five are listed 



in the compilation by Kroll & Kock (1987), and all five 
are found in the work by Moity (1983). Unfortunately, 
there is a systematic difference of ~0.20 dex between these 
two investigations (Moity's values being lower). This has 
some consequences on the final iron abundance derived 
from these (few) lines, and certainly has a non-negligible 
effect on the derivation of log g from the ionization bal- 
ance. We ran a few tests and reached the ionization bal- 
ance for log ff=4.3 if mean ^/-values were adopted, and 
for log g=4.1 if Moity's values were used. These results, 
together with the paucity of Fe n lines, made us decide to 
keep log g= 4.1 ± 0.2 as our final choice for the analysis. 

4. Abundance Analysis Highlights 

The analysis of the beryllium spectral region was per- 
formed via spectrum synthesis. For this purpose, we 
adopted a model atmosphere computed with T c g — 
6400 K, log g= 4.1, and [Fe/H] = -3.3. Microturbulence 
was assumed to be 1.5 km s . 
The ID spectrum synthesis analysis was carried out with 
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Fig. 1. Reduced and normalized observed spectrum of 
G 64-12 with overplotted our best-fit synthesis (thin line) 



Kurucz model atmospheres and SYNTHE (Kurucz 1993|) . 
The model was convolved with a Gaussian function of 
FWHM = 6.7 km s _1 (which combines instrument, ro- 
tation and macroturbulence broadenings), and was com- 
puted with the a-elements enhanced by +0.4 dex and 
with the "approximate overshooting" option switched off. 
The solar abundances were taken from the compilation of 
Grevesse & Sauval (1998). The list of lines necessary to 
synthesize the near-UV spectral region is the one tested 
by Prirnas et al. ( |1997| ). 

Additionally a time-dependent 3D hydrodynamical model 
atmosphere with parameters appropriate for G 64-12 
has been constructed, with the same 3D, compressible, 
radiative-hydrodynamics code which has previously been 
successfully used for studies of solar and stellar convection 
(e.g. Stein & Nordlund |1998| Asplund et al. [1999| , |2000| ). 
The equations of mass, momentum and energy coupled 
to the 3D radiative transfer equation (which includes the 
effects of line-blanketing) have been solved on a Eulerian 
mesh with 100 x 100 x 82 zones covering 21 x 21 x 8.5 Mm. 
For the 3D LTE spectrum synthesis, a representative se- 
quence of 1 hr stellar time with snapshots every 1 min was 
selected; the temporal average of T g for the sequence is 
6460 ± 27 K, i.e. very close to the value indicated by the 
IRFM method. It is noteworthy that no free parameters 
like mixing length parameters, micro- and macroturbu- 
lence enter the construction of 3D model atmospheres or 
spectral synthesis, as the time-dependent convective mo- 
tions and inhomogeneities are self-consistently calculated. 
In order to quantify the systematic errors introduced when 
relying on ID analyses, a strictly differential comparison 
of the 3D predictions have been carried out with a corre- 
sponding MARCS model (Gustafsson et al. 1975) with the 
same parameters and input data (Asplund et al. 1997). 



4-.1. On other Elemental Abundances 

Because of the large spectral coverage (see Table 1), all the 
most commonly studied elements (a, iron-group, heavy) 
are accessible from the observed spectrum of G 64-12. Al- 
though the complete abundance analysis of G 64-12 is be- 
yond the scope of this work, some of the elements most 
relevant to it deserve mentioning. 

a-elements: we checked the abundance of magnesium and 
calcium in order to test how representative G 64-12 is of 
the population of the very metal-deficient stars. From the 
equivalent widths measured for two Mgib lines and four 
Cai lines we derive [Mg/Fe] = +0.43 ± 0.03 dex and 
[Ca/Fe] = +0.48 ± 0.07 dex, i.e. "normal" and in good 
agreement with what can be found in the literature. 
Oxygen (from OH lines): the same spectral region where 
the Ben doublet falls is also very rich in OH molecular 
lines, which represent one of the few oxygen abundance 
indicators available to spectroscopists. Our ID LTE spec- 
trum synthesis analysis suggests [O/Fcl = +1.1, similarly 
to wha t foun d by Israelian et al. ([1998) and Boesgaard 
et al. ( 1999a ) using the same abundance indicator and 
for stars of comparable metallicity. Our finding is also in 
agreement with the recent oxygen determination made by 
Israelian et al. (200C) in the same star. However, due to 
severe systematic errors affecting the near-UV OH lines 
in metal-poor stars, we prefer not to adopt this value. We 
suspect that the oxygen abundance derived from the UV 
OH lines is in error, being affected by metallicity depen- 
dent 3D and/or NLTE effects (Asplund et al. [1999D - A 
future analysis of Be in a larger sample of stars will deal 
with this issue (Prirnas et al., in preparation). 
Lithium: the S/N ratio achieved around the Li I line at 
670.8 nm is « 300. Our measurement of the Lil 670. 8nm 
equivalent width gives EW = 22.2 ± 0.5 mA, in agreement 
with the value published by Ryan et al. ( |1999| , EW = 21.2 
± 1.1 mA). The lithium abundance derived is A(Li) = 
2.23 ± 0.01 dex. The difference with Ryan et al. fll999| ) 
value (A(Li) = 2.14 ± 0.03 dex) can be explained by the 
difference in adopted temperatures. 

5. The Beryllium Abundance and Its Uncertainty 

The best-fit synthesis of the Be n doublet is achieved with 
log N(Be/H)= -13.20 for the 313.0 nm line, and with 
log N(Be/H) = -13.10 for the 313.1 nm line. Being both 
lines well resolved, equal weights were assigned. The beryl- 
lium abundance thus derived using Kurucz ID LTE model 
atmosphere is log N(Be/H)= —13.15 ± 0.15 dex, signifi- 
cantly higher than what is expected from extrapolating 
the previously observed trends to these low metallicities 
(under the assumption that the Be abundance keeps de- 
creasing) . The magnitude of the error-bar assigned to our 
measurement mainly comes from the dependence of the 
beryllium abundance on changes in the stellar parameters. 
In order of importance, we find that ABe = ±0.1 dex for 
Alog g = ± 0.2, ABe = ±0.04 dex for AT cff = ± 100 K, 
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Fig. 2. An enlarged view of the Bell resonance doublet in 
G 64-12 (histogram), together with our best-fit synthesis 
(thin line). Overplotted are two syntheses computed with 
the beryllium abundance increased and reduced respec- 
tively by 0.15 dex (dotted lines) 



and ABe = ± 0.03 dex for A[Fe/H] = ± 0.1 dex. To 
these values (to be summed in quadrature) one has to add 
the dependence of Be on the placement of the continuum. 
This uncertainty in such metal-poor star is estimated to 
be ±1% at most, implying a variation in Be of ±0.05 dex. 
Fig. 2 shows an enlarged view of the Be doublet, with our 
ID LTE best-fit synthesis, and two others computed with 
± 0.15 dex variation in the Be content, which demonstrate 
our detection and measurement accuracy. 

As expected, the predicted Ben line strengths with the 
3D model atmosphere are quite similar to the ID case, be- 
cause the lines are formed in the deep atmospheric layers 
due to the lines being very weak and coming from the 
majority ionized species. Assuming LTE, the abundance 
correction derived by using 3D hydrodynamical model at- 
mospheres is only 0.03 dex for Ben 313.1 in G 64-12, in 
the sense that the ID abundance should be corrected up- 
wards with 0.03 dex. It should be emphasized that the 1D- 
3D comparison has been done strictly differentially. The 
quoted differences are therefore estimates of the system- 
atic errors due to the adopted model atmospheres one can 
expect when using a ID analysis, regardless of whether 
Kurucz or MARCS model atmospheres have been used. 
Although the above-mentioned 3D effects are small, it 
should be borne in mind that these estimates assume LTE 
and there is no guarantee that the same will be true for 
NLTE. Recent improved ID NLTE calculations appropri- 
ately done for G 64-12 with a very extended Be atom (71 
levels compared to the 9 used by Garcia Lopez et al. 1995| ), 
find some small but significant NLTE effects. In terms of 
abundance, the NLTE abundance is 0.07 dex larger than 
the LTE result. In spite of the similarity with the LTE 
result, it should be noted that the line formation for the 
Bell lines is far from being in LTE, since the dominant 



NLTE effects - over-ionization and over-excitation - have 
opposite influence on the line strength and thus the effects 
partly cancel. This may suggest that for some other stellar 
parameters the NLTE effects could be significantly larger. 
The details of the ID NLTE calculations for Be n appro- 
priate for a range of stellar parameters will be presented 
elsewhere (Garcia Perez et al., in preparation). It also re- 
mains to be investigated whether the steep temperature 
gradients and inhomogeneities in 3D model atmospheres 
may drive more pronounced departures from LTE, similar 
to re cent calculations for Li and O (Asplund & Carlsson 
[2000; Asplund et al., in preparation). 



6. The [Be/Fe] ratio 

Our claim to have detected a significantly higher beryllium 
abundance in G 64-12 (compared to other low metallicity 
stars) can be more clearly seen when our new (ID LTE, 
for consistency with the rest of the sample) determination 
is plotted in a log N(Be/H) vs. [Fe/H] graph (cf Fig. 3a) 
together with carefully selected literature data of simi- 
lar quality. All the data shown in Fig. 3 (both panels) 
had their Be abundances computed assuming a "high" 
T e g scale, hence consistent with what we did. 



Extrapolating the Boesgaard et al. (1999b) Be vs Fe 
relation at lower metallicities, our ID LTE Be value is 
0.6 dex higher (using their single line fit) or 0.8 dex higher 
(using their double line fit), regardless of the T e g scale 
adopted. One may also compare G 64-12 to BD-04 °3208 
(similar T c g and log g but a factor of 10 higher metal- 
licity, cf Boesgaard et al. 1999b) and note that their Be 



contents differ by a factor of 3 only (instead of the ex- 
pected factor of 10, should Be and Fe grow at the same 
rate). From our measurements we also infer [Be/Fe] = 
2.15 ± 0.19, which is at least 0.5 dex off the average 
trend (Fig. 3b). Although we must say that G 64-12, to- 
gether with the recent Be detection in LP 815-43 ( [Fe/H ] 
= —2.95, open triangle in Fig. 3; cf Primas et al. |2000| ), 
gives more weight to the hypothesis of a possible flattening 
of the Be trend at the lowest metallicities, at the moment, 
the interpretation of our new result is not straightforward. 
According to Vangioni-Flam et al. (1995), the finding of 
high Be abundances in stars with metallicity below [Fe/H] 
= — 3.0 may favor shock acceleration in the gaseous phase 
of superbubbles produced by collective SNII explosions as 
the main mechanism of Be production. This, in turn, im- 
plies that the major role is played by the most massive 
stars only (say, with initial mass M > 60Mq). 
According to the Parizot & Drury ( 1999] ) scenario, a 



double trend of Be vs O (and Fe) is expected depend- 
ing on if the star formed inside a superbubble (hence 
with a high Be content due to the effect of collective 
SN explosions) or outside. In this case, we cannot ex- 
clude the possibility that G 64-12 may be representa- 
tive of the "high" trend, whereas other previously studied 
very metal-deficient stars (like BD— 13°3442, for instance, 
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Fig. 3. Panel (a): log N(Be/H) vs. [Fe/H] for the high- 
est quality data points currently available, with the Boes- 
gaard et al. ( 1999b| ) least squa re fits overplotted (open 
squares from Boesgaard et al. ( 1999b), th e open trian- 
gle is LP 815-43 from Primas et al. ( 200C ), whereas the 



filled symbol represents this work - our ID LTE value, 
for consistency with the rest of the data points); Panel 
(b): [Be/Fe] vs [Fe/H] for all the stars plotted in Panel (a) 
(same symbols have been adopted) 



Boesgaard et al. 1999b ) formed out of gas outside a super- 
bubble. In order to test this hypothesis, new detections of 
beryllium in this challenging metallicity range are needed. 
One last possibility that should not be underestimated is 
the effect due to 3D and NLTE corrections. As previously 
stated, NLTE corrections have usually been assumed to 
be negligible, but they could be important for different 
combinations of stellar parameters. The future results of 
Garcia Perez et al. work will shed light on this issue. 

7. Conclusions 

We have presented a new analysis of the near-UV part 
of the spectrum of G 64-12 with the aim of studying its 



beryllium content. The high quality of this UVES (Sci- 
ence Verification) spectrum has led to a clear detection 
of beryllium. This is the first measurement of Be in a 
star of such low metal content ([Fe/H] = —3.3). We have 
performed our analysis with ID LTE model atmospheres, 
but introduced the novelty of correcting it for both 3D and 
non-LTE effects. Hence, our final 3D NLTE Be abundance 
is log N(Be/H) = -13.05±0.15. We have shown that when 
compared in a consistent manner to previous high quality 
literature data, the Be abundance here detected is signifi- 
cantly higher than what expected. The strength of the OH 
molecular lines is comparable to what observed in other 
stars of similar metallicity, thus excluding the possibility 
that the high Be content measured in G 64-12 may be 
somehow related to an abnormally high O content. Fur- 
thermore, our measurement disagrees with the statement 
of a parallel (i.e. same rate) increase of Be and Fe, which, 
in turn, suggests a flattening or the presence of dispersion 
at given metallicity in the Be vs Fe trend during the early 
evolutionary phases of our Galaxy. 
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